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Abstract 

The development and application of a multi-dimensional capability for modeling and 
simulation of aviation-sourced particle emissions and their precursors are elucidated. 
Current focus is on the role of the flow and thermal environments. The cases 
investigated include a film cooled turbine blade, the first-stage of a high-pressure 
turbine, the sampling probes, the sampling lines, and a pressure reduction chamber. 


1. Introduction 

Aircraft emissions can influence the climate change of the atmosphere and health 
responses in the vicinity of the airports. To alleviate these increasingly significant 
impacts on the environment, it is necessary to acquire a better understanding of the 
formation and subsequent development of gaseous pollutants, aerosols (volatile 
particles and soot) and their precursors in the internal flow and in the plume of the jet 
engines operating over the full range from ground to flight altitude. This could be 
achieved through efforts which combine the use of accurate measurement and 
high-fidelity modeling and simulation. 

The volatile nature of the aerosols introduces difficulties in their repeatable and 
accurate measurement because their concentration and size can be quite sensitive to the 
sampling conditions and procedures. These difficulties are sometimes further 
compounded by the necessity of in-situ measurements, i.e., sampling from aircraft in 
flight. Clearly, the ability to predict the effects of sampling conditions and procedures 
on the aerosol behavior will greatly assist the development of sampling techniques and 
the optimization of sampling procedures for making accurate and quantitative 
measurements of particulate emissions. 

Current modeling studies often have used low-order methods, where chemical kinetics 
and/or particle microphysics are driven by averaged flow parameters specified as a 
function of time (e.g. [1], [2], [3], [4]). In some flow regions and for certain chemical or 
microphysical processes, these low-order methods could adequately capture the overall 
activities. Nevertheless, the more complex multi-dimensional analysis needs be applied 
to flow regions in which the chemical or the microphysical activities could be 
appreciably influenced by flow non-uniformity, particularly in thermodynamic 
conditions and residence time. 
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This paper describes a Computational Fluid Dynamics (CFD) based modeling and 
simulation tool for complementing limited measurement studies and addressing the 
sampling system issues. It accounts for and integrates fluid dynamics, chemistry and 
particle microphysics relevant to aircraft emissions. The discussion proceeds as follows. 
Section 2 presents the conservation equations for the gas-phase flow. Section 3 
describes the continuity equation for aerosol transport. Section 4 discusses the 
implemented models of various microphysical processes such as thermophoresis, 
coagulation, nucleation, soot activation, condensation, and the interaction between the 
gaseous species and the particles. Section 5 presents the application of this 
multi-dimensional modeling and simulation tool. The cases investigated include a film 
cooled turbine blade, the first-stage of a high-pressure turbine, the sampling probes, the 
sampling lines, and a pressure reduction chamber. Section 6 summarizes the current 
state of the progress and suggests areas warranting future improvements. 


2. Governing Equations for Gas-phase Flow 

The conservation equations for the gas-phase flow in the multiple rotating frames of 
reference (MRF) which is convenient for including the rotational components of the jet 
engine in the computational domain are summarized in a hybrid formulation below. 
Continuity equation: 


Species transport equation: 

fill A dv * 1J>> -o,y^ ■ dA + JJJ fjsv 


( 2 ) 


Momentum equation: 


d_ 

dt 


JJJ pudV- 1- JJ pu(u-U g )dA--jj pdA + JJ r dA 


-JJJ/* 


xudV 


(3) 


Energy equation: 


EdV + \\ pE(u -U g )-dA = -JJ pu ■ dA - JJ ? ■ 

rr 


dA 


(4) 


Turbulence equations: 
d_ 
dt 


J t \\. f dV + \\ A pk{u -U g )-dA = S k + \\y + ■ dA 


(5) 
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d_ 

dt 


\\\vP £ dV + \\ g ) dA = S £ + ^ (p + —)Ve dk 


( 6 ) 


where 

q = -ldVT, 


f = ~/i(V-u) I + ^-//(Vu + (Vu) T ), 

V is an arbitrary control volume with control surface A, p is, the fluid density, p m is 
the partial density of species m , w m is the production rate of species m per unit 
volume, u is the flow velocity in a stationary Cartesian coordinate system, a), is the 
grid rotating velocity for each designated axis, U =© s xr, E is the total energy 
density, k is the turbulent kinetic energy, e is the turbulent dissipation rate, S k is the 
source term for k , S e is the source term lor e , p is the static pressure, T is the 
viscous stress tensor, q is the heat flux vector, D tJf is the effective mass diffusivity 

coefficient which is equal to , k , p and p, are the thermal conductivity, laminar 

viscosity and turbulent viscosity, respectively. For gas, the Schmidt number ( Sc ) is 
about 1. Different angular velocities (and even different rotating axes) are assigned to 
different mesh blocks or groups within the model. Balance equations for each group are 
expressed in the relative reference frame but in terms of the absolute velocity (i.e. the 
velocity with respect to a stationary coordinate system). 


3. General Dynamic Equation 


A general dynamic equation (GDE) is the continuity equation for aerosol transport. 
Because the GDE is a nonlinear, partial integro -differential equation, analytical 
solutions are available for only a few special cases [5, 6]. Aerosols are often described 
in terms of size distributions, i.e. the number of particles in a given size bin. Therefore a 
multiple size bin model is adopted here to approximate the aerosol general dynamic 
equations (GDE). Since the governing equations for gas-phase flow are written in the 
multiple rotating frames of reference (MRF), the GDE will be cast in this reference 
system too. Furthermore, it is written in a discrete form as a population balance for each 
cluster or particle size, and it describes particle dynamics under the influence of various 
physical phenomena: convection, diffusion, coagulation, nucleation, condensation, 
evaporation and external force field. 


d_ 

dt 


JJJ A M M = \\(D„ 


+ 


2v, 

3 VcZ 


)VN k -dA 



dN_ 
v dt 


k 


dV, 

micro 


k = l,---,M B 


(7) 


where N k is the number density of particles in size bin k , M B is the total number of 
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the size bins, u is the flow velocity in a stationary Cartesian coordinate system, 1 v, is 
turbulent kinematic viscosity, C u is 0.09 , f u is the damping function of the turbulence 
model. The effective molecular (Brownian) diffusivity coefficient, D M , is defined as 


D, 


k B TCC k 

3 m g d pjc ’ 


( 8 ) 


All the variables in the above equation will be defined in the following sections, 
represents the summation of microphysical processes, 

I micro 


dN± 

dt 


•V 

II 

+ dN k 

dN k 
+ — L 

+ dN k 

dt \ 

\ micro 1 

drifting 4t 

coagulation 

nucleation 


(9) 


The first term in Equation (9) is associated with the drifting velocity caused by external 
forces, e.g., the thermophoretic force, the lift-drag force, the electrophoresis and 
gravitational force. As for the rest of the terms, coagulation means that the particles 
collide and attach to each other; nucleation represents the formation of the nuclei of 
new particles; and condensation indicates that molecular mass transfer from the vapor 
phase to the particles, while evaporation is the reverse process of the condensation. In 
the present work, the particle number density on a solid wall is assumed to be zero or a 
fraction of the value at a point next to a solid wall. 

4. Models for Particle Microphysical Processes 

The most important and efficient binary nucleation in the atmosphere is that of sulfuric 
acid and water. A free sulfuric acid molecule tends to gather water molecules around it 
to form hydrates. The main reason is that H2SO4 has extremely low saturation vapor 
pressure. Doyle [7] showed that even when relative humidity is less than 100% , 
extremely small amount of H 2 SO 4 vapor are able to induce nucleation. Within the 
turbine engine exhaust, major gas-phase chemical reactions will take place, and soot 
activation [8, 9] into water condensation nuclei also occurs. If the surface of a foreign 
particle (e.g., soot) does participate in the gas-to-particle conversion process, the 
nucleation is referred to as heterogeneous. More detailed discussions on soot activation 
and condensation of water and sulfuric acid can be found in Wang’s work [9]. The first 
elemental step in heterogeneous nucleation involves molecular adsorption which is 
between soot and H 2 SO 4 and SO 3 . The nano-size aerosols will collide with each other 
and stick together. Such a process is the called coagulation and it will change the 
particle size distributions. 

4.1. Characteristics of Particles 

In the present work, we assume that 

• All particles are assumed to be spherical. 

• Heat and momentum transfer between the carrier gas and the particles is 
negligible. 

• Gas turbulence influences the particles, but the particles do not 
influence the gas turbulence. 


'For noninertial particles (i.e. extremely small aerosols) advected by a turbulent flow, 
particle velocity coincides with carrier velocity. In the present work, all particles are 
treated as noninertial for simplicity. 
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No particle is charged. 

The particle size distribution of soot entering the post-combustor region 
is always prescribed. 

All particles are noninertial. 


Based upon the sectional approach, the total volume range for the particles is divided 
into M b bins. Using a geometric progressing factor, the volume of a particle in each 
bin can be calculated easily. For example, in an aerosol system, the particle size ranges 
from 1 run to about 1 pm, hence, the corresponding volume ranges from 10 21 nr to 
10“ l8 m 3 . To cover the 9 orders of magnitude volume span with 12 bins, the bins are 
linearly spaced on a logarithmic scale, so that v k = l0 9l(12 ~ i) v k _ l , where v k is the 
volume of a particle in size bin k. 


4.2. Thermophoresis 


Thermophoresis arises from the temperature gradients in the carrier gas: small particles 
drift from hot regions towards cold regions. Brock [10] derived the following 
expression for the thermophoretic velocity, y TH 


V TH, 


-3Cc k (k„ + c,k p Kn k ) 


0 + 3 c m Kn k )(k+2k + 2c t k v Kn k ) 


}h_ 


\VT 


( 10 ) 


where //„ is the viscosity of the carrier gas, k g and k p are the gas and particle thermal 


conductivities, respectively, and c m and c t are phenomenological coefficients known 
as the ’’isothermal slip” coefficient and the ’’temperature jump” coefficient, 
representative values according to Brock [10] are c t - 2.2 and c m = 1.0 .The smaller 

the diameter of a particle, the larger the effect of the thermophoretic force. Since the 
turbulent thermal diffusion discussed by Elperin [11] is not included in the present 
model, we have 


dt 


= - V -(VmA) 


I drifting 


(ii) 


When the Knudsen number is much larger than one, i.e., Kn k » 1 (particles are much 
smaller than the mean free path of the carrier gas), an alternative form for the laminar 
thermophoretic velocity is given by [12]: 


V TH,. 


4{\ + 7rod%) 


PJ 


\VT 


(Kn k » 1) 


( 12 ) 


where the negative sign indicates that the motion is in the direction of decreasing 
temperature, and a is the nondimensional thermal accommodation coefficient, its 
typical value is about 0.9. The thermophoretic velocity for Kn k » 1 is independent of 
the particle size. Furthermore, it is nearly independent of the particle material. 

4.3. Coagulation 

Aerosol coagulation is important because it changes the size distribution and the 
composition of particles. Coagulation is a process where particles in a population 
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collide with each other and stick together to form larger particles [13]. The total volume 
is conserved in a coagulating aerosol ensemble. Collisions between aerosols can be 
initiated by gravitational settling, turbulence, and thermal (Brownian) motion. For 
particles smaller than one micrometer in diameter, the most important collision process 
is the Brownian motion, which is considered in the current work. Different approaches 
have been developed in the past to simulate coagulation, depending upon the need and 
availability of computer resources. Since the model has to be incorporated into a 
three-dimensional CFD program, a semi-implicit scheme [14] is adopted here. All the 
particles are assumed to be spherical. Soot density is less than that of carbon black and 
usually in the range of 1700-1800 kglm 3 depending on the porosity of soot. Soot 
particles are generally small, ranging in size from 5 nm to 90nm, but may be up to 
several micrometers in extreme cases. Soot scavenging by coagulation occurs when a 
small H2SO4-H2O droplet collides with a larger coated soot. 


4.3.1. Single particle size distribution 


The rate of change of a single particle size distribution (PSD) due to coagulation is 
given by a modified Smoluchowski Equation, 


dN, 

dt 


=7EVA-/,r»T s .A 

' Z 2=1 2=1 


( 13 ) 


where B k . ( m 'Ks particle) ) is the coagulation kernel or collision frequency function 
of two colliding particles. A volume-conserving solution [15] is given by 


v,. 


v k N kJ = ■ 


n >, a 




( 14 ) 


where v k is the volume of a particle in size bin k , M B is the total number of bins, h is 
the CFD time step used in the solution of the carrier gas. When two particles collide and 
stick together, the volume of the intermediate particle is 

V iJ = V i +V j ( 15 ) 

If this volume falls between two bins, the new particle is split into adjacent bins as 
follows. 


fi,i,k = ■ 


V * +I ~ v u 

_v k+l -v k 

1 

0 


v >.j 


V k ^V i , j <V k+ 1 k<M B 

Vk-i<Vi.j<Vk k>] 

Vi,j>V k k=M B 

all other cases 


( 16 ) 


The generalized Brownian coagulation kernel B i . for collision of particles of volume 
v i and Vj is given by Fuchs [13] 


2 ^d i +d J )(D pi + D pJ ) 

dj+dj 8< £>,,+£>,,,) 

d +d j +(4 2 +^j )" 2 ( d +dj m cl,t+c 2 p.jf 12 


( 17 ) 
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where d i and d ] are the diameters of particles in the i th and j th size bins, the unit of 
B j j is mils particle ) , D p j or D p ; is the particle diffusion coefficient in the i th 
or j th size bins and is given by 


k B TCc i 

3 xd pJ n g 


(18) 


where k B is the Boltzmann constant (= 1 .380658 x 1 0 2 : JK 'molecd' ), T is the 
absolute temperature (K), d pj is the diameter of the size-i particle, and ju is the 
dynamic viscosity of the carrier gas. The Cunningham-slip factor is defined as 

Cc, = 1 + Kn t (\ .257 + 0Ae ' ' IKn ‘ ) (19) 


where Kn i is the Knudsen number defined as the ratio of the mean free path of the 
carrier gas, k , to the diameter of the size-i particle 

2k 


Kir = 


d. 


(20) 


The mean free path of the carrier gas depends upon the pressure, temperature, 
molecular weight of the carrier gas, and the universal gas constant, it is given by 


\=- 


P 


( 8 wd ' 

A ) 


(21) 


„ j is the mean thermal speed of a particle in size bin i, 


C p,i 


/ \ 1/2 

^ 8 k B T ' 

V Km < J 


( 22 ) 


where m i is the mass of a particle in size bin i . S ! is a correction factor which 
represents the mean distance between the center of a sphere and a particle which has 
bounced from the surface of the sphere and has traveled a distance of one particle mean 
free path, it is given by 


(d rj+ k pj y-(d; j+ ^r 

3d. /. 


The particle mean free path of a particle in size bin i is 

8 D 


' Kc P ,i 

Finally, the change of the size distribution due to coagulation is 


dN, 

dt 


N ka ~N kJ _ h 


(23) 


(24) 


(25) 


4.3.2. Multiple particle size distributions 

Coagulation among multiple particle distributions (e.g. among the soot particles and the 
H 2 SO 4 -H 2 O liquid droplets) is also considered in the current study. Coagulation among 
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two independent size distributions will add one more particle size distribution for the 
resultant particle, hence, the total number of particle size distributions becomes three. 
Here, the semi-implicit coagulation solution is extended to treat coagulation involving 
three different size distributions (e.g. N T = 3), with size bins in each distribution being 


M b . The resultant particle number density of the Y distribution in bin k at time t 
during coagulation is defined as 


N 


Yk.t 


N Yk,t-h + h ( T 1 +^>) 

1 + hT 3 


(26) 


where 




V, 


Yk M= 1 


P y,mYj 


N Mj.t-h fYi,Mj,Yk P Yi,Mj^Yi^Yi,t 


(27) 


‘ 1 1 




Yk M = 1 /= 1 


Ql M K / , 
7=1 


N Mj.t-h f Ii.Mj. Yk P Ii.Mj^Ii^Ii,t 


(28) 


m b 

^ =z 


7=1 


N t 

Z[(l -b,M XI fYk,Mj.rk^ + P Y,M^ P Yk,Mj^Mj,t-h > 
M= 1 


(29) 


Term T x accounts for the production of larger particles in distribution Y from 
self-coagulation and from heterocoagulation between distribution Y and distribution 
M . Term T 2 accounts for the production of particles in distribution Y from 
heterocoagulation among the other two independent distributions ( I it Y and M &Y). 
The first part of term T, accounts for self-coagulation loss in distribution Y to larger 
sizes, while the second part accounts for loss of distribution Y to all other distributions 
due to heterocoagulation between the Y distribution and the other distributions. 


If the volume of the new particle falls between two bins, it is split into adjacent bins as 
follows. 


f, 


Ii,Mj,Yk 



v -V 

u Yk + 1 v Ii,Mj 

v Yk 

P Yk ~ P Ii.Mj < Y>rk + 1 

k<M B 



V 

V Ii.Mj 


_ V YM~ V Yk 


{ 1- f 

A J li.Mj.Yk- 

-1 


k> 1 

1 


P Ii.Mj — 

k=M B 

0 


all other cases 



(30) 


The values of P , Q , and L in equations (27-29) are either 1 or 0, depending on the 
specifics of the coagulation considered. Parameter P Y M = 1 , if the coagulation between 
particles in distribution Y and in distribution M produces larger particles in 
distribution Y . For example, P soot soot = 1 , and P HlSOi -H 1 o,H 1 so t -H 1 o = 1 • Parameter 
Qi my = 1 , if the coagulation between particles in distribution I and in distribution M 
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produces particles in distribution Y , where I ^ M and I i^Y . For example, 
Qsoo,,H 2 so i -H 2 o,mi x = 1 > and Q HiSOi _ HiOMix = 1 • Parameter L, M = 1 , if the coagulation 
between particles in distribution Y and in distribution M does not produce particles in 
distribution Y . For example, L oot H ^_ Hp = 1, and L HiS0i _ Hi0 soot = 1 . 


Finally, the change of the size distribution for each type of aerosol due to coagulation is 


dN v 


dt 


N - N 

ly Yk.t ly Yk,t-h 


(31) 


4 . 4 . Homogeneous Binary Nucleation of H2SO4-H2O 


The classical nucleation theory assumes that as the state of current phase is becoming 
unstable, a new phase is formed. It is based on the equilibrium state of small liquid 
drops, or embryos, in contact with their parent vapors; some of the embryos, which 
become freely growing droplets, are said to be "nucleated". The theory predicts the 
nucleation rate as follows: 


where 4 hom 


4 hom ~?td* Ph.SO. ^ ICO exp 


f ag* a 


V RJj 


is the homogeneous nucleation rate ( particle! (m 3 s ) ), ltd* 2 fin, so 0 


(32) 
is a 


prefactor term, AG* is the change of the Gibbs free energy ( J /mol ) during the phase 
change, R u is the universal gas constant. The change of the Gibbs free energy is a 
function of the number of sulfuric acid and water molecules and is represented by a two 
dimensional energy surface. To become stable, a cluster will have to follow the path of 
least energy on this surface, thus leading to the so-called saddle point. To avoid solving 
iteratively to find the roots of the equations which satisfy the saddle point condition, the 
parameterization of Kulmala et al. [16], updated by Vehkamaki et al. [17] for low 
temperature emissions and by Vehkamaki et al. [18] for high temperature emissions, is 
used here (See Appendix A). This scheme provides the nucleation rate J H S0 __ H 0 

(number of new embryos formed per second and per cubic meter) in the water/sulfuric 
acid mixture, and the critical cluster composition (total number of molecules and 
sulfuric acid mole fraction). These parameterized equations reduce the computing time 
by a factor of 500 compared to non-parameterized nucleation rate calculations. 


4.4.1. Low Temperature Range 

The parameterized formulas are valid for temperatures between 1 90. 1 5 K and 305 . 1 5 
K, relative humidity between 0.01% and 100%, and total sulfuric acid molecule 
number from 10 10 to 10 17 in 3 . The parameterization is limited to cases where 
nucleation rates are between 10 1 and 10 16 in V and the critical cluster contains at 
least four molecules. For the cases where the temperature is below 190. 15 K , a zeroth 
order extrapolation is applied to the nucleation rate. More details can be found in 
Appendix A. 
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4.4.2. High Temperature Range 


The parameterized formulas are valid for temperature between 300.15 K and 400.15 
K , relative humidity between 1% and 100%, and total sulfuric acid molecule number 
from 2xl0 15 to 5xl0 21 molecule m . The parameterization is limited to cases where 
nucleation rates are between 10 5 and 10 2 ° particle m 'A and the critical cluster 
contains at least four molecules. For the cases where the temperature is above 400.15 
K , a zeroth order extrapolation is applied to the nucleation rate. More details can be 
found in Appendix A. It is noted here that, for sulfuric acid solution, experimental data 
are available for its surface tension when T <323 K , for its density when T <313K , 
and for its activities 2 when T <350 K . The computed nucleation rates and other 
thermodynamic properties may be of greater uncertainty when the temperature is out of 
the experimental range. 

4.5. Soot Activation 

Since soot is hydrophobic, it must be activated in order to be capable of taking up water 
from the gas phase. Soot can be activated as water condensation nuclei by adsorption 3 
of oxidized sulfur species (H2SO4 and SO3) and by scavenging, i.e., coagulation with 
volatile sulfate aerosols. 

The fraction of soot surface due to the above two activation pathways is 0 ! = 0 ads , + 
0 sca ( , where i is the index of the size bin. Karcher [8] presented a kinetic formula for 
0 

w ads,i * 

^±£ = .25 ac(C, 0 +C H SY , )x 6.02x10 23 (33) 

dt 503 lSOi <7 0 


where a is the sticking probability. It is either set to zero for T > 420 K or set to 
unity for T < 420 K . c is the mean thermal speed for SO3, C is the molar density of 
gas phase species, ( mol tin ), and <r 0 is the average number of sites per unit area of soot 
surface (< 7 0 = 5x1 0 18 m ~ 2 ). This formula which assumes that the vapor adsorption on 
fresh soot takes place in the gas kinetic regime gives the maximum adsorbed sulfur 
mass, and it represents the upper bound of the binary H2SO4-H2O heterogeneous 
nucleation rate on the soot surface. 

The change of surface coverage for a soot particle in size bin i by scavenging volatile 
(H2SO4-H2O) droplets from size bin j = 1 to j — M n is determined by 


dO ■ £4 , 1-0 

— S -^ L = 25n x V fi. NA- x — r— L 

dt fr 7 7 7 7 << 


( 34 ) 


where B i . is the coagulation kernel between soot in size bin i and H2SO4-H2O in size 


2 The product of the mole fraction of a solution component and its activity coefficient is 
defined as the activity of the component. 

3 When a gas is brought into contact with a solid particle and part of it is taken up by the 
solid surface, the process is known as adsorption. If the interactions between the gas 
and the solid are weak, such as between sulfur species and soot, the process is called 
physical adsorption which is characterized by the molecular bonding. 
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bin j, N j is the number density of H 2 SO 4 -H 2 O in size bin j, d 2 . is the diameter of 
H2SO4-H2O droplets in size bin j , and d 2 , is the diameter of soot in size bin i . 


The gas phase depletion of H2SO4 and SO3 due to adsorption is proportional to ( 1 - ), 
i.e., the available "dry" soot surface area, they are determined from as follows. 


dC 


so. 


dt 


= -.25acC SO) Y J N i d 2 pJ (\-0 i ) 


(35) 


dC„ 


dt 


= -.25acC HiSO £N i d 2 p . (l-Of 


(36) 


where Af is the number density of soot in size bin i. This model also accounts for the 
condensation of H 2 SO 4 -H 2 O on the "wetted" soot surface area. Consequently, 
condensation rates onto soot particles and corresponding gas phase depletion of H2SO4 
and H 2 O are proportional to 0 I . 

Finally, using a semi-implicit method, Q i is evaluated by 

c ,_ er h +h(Q ads +Q sca ) 

1 + HQ ads + Q sca ) 

where 

Q ads = .25 — ci (C, 0j + C HiSOi ) x 6.02 x 10 23 (38) 

<y n 


and 


Q =.25^xY B N.d 2 . 

SCsca j2 l.J J P.J 

a p.i 7=1 


(39) 


4.6. Condensation of Water and Sulfuric Acid 


If the vapor pressure of H2SO4 or FFO is in excess of the equilibrium vapor pressure, 
condensation will occur. Condensation can occur on the liquid H2SO4-H2O aerosols as 
well as on the liquid coating of soot surface. If the vapor pressure of H2SO4 is less than 
the equilibrium vapor pressure, then the sulfuric acid molecules will evaporate from the 
H2SO4-H2O aerosols. However, the equilibrium vapor pressure of H2SO4 over the 
liquid droplets is so small that, once H2SO4 molecules are condensed, they can hardly 
re-evaporate. In contrast, H 2 O will condense to and evaporate from the droplets 
depending on the evolution of temperature and the partial pressure of H 2 O in the carrier 
gas. 

Following Fukuta and Walter [9, 19], the condensational growth of m-H 2 0, or 


H,S0 4 onto particles in the size bin k is given by: 


dn„ 


dt 


■ ~ 4ftr k D* vjn . k 


Pm P m.k.sat 

RT 


(40) 
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where is the rate of change in moles of species m due to condensation onto or 
evaporation from an aerosol with radius r k , the modified diffusivity is 


D 


I^vjn 

r, 4 D r 


(41) 


where X is defined in Equation (21), D v m is the diffusivity of species m, as defined in 
Appendix B. 

The mean thermal velocity of gaseous species m is given by 


r ni /2 

' 8 RT ' 


C m 


7TMW 


(42) 


m J 


It should be noted here that, when calculating the condensation on soot particles, since 
only the activated surface of soot can accept the diffusion of the molecules of H 2 SO 4 
and H 2 O, the right hand side of equation (40) needs to be multiplied by the fraction 9 k 


provided via equation (37). 

The saturation vapor pressure is computed from Kelvin equation which implies that the 
curvature will increase the value of the saturation vapor pressure of species m, 


Pm,k,sat 

n flat 

r m,sat 


- exp 


2<7v 


V R u Tr k J 


(43) 


where <r is the surface tension of a liquid droplet ( Joulelm 2 ), y m is the partial molar 
volume of species m in a solution ( in 3 / mol ), p m k sat is the saturation vapor (partial) 
pressure of gaseous species m over a spherical droplet of radius r k , and p; n n a, sat is the 

saturation vapor (partial) pressure of species m over a flat liquid surface with the same 
composition as the droplet. The effect of curvature is important for droplets smaller 
than 0. 1 pm . Further details of various thermodynamic properties of aerosols are given 
in Appendix B. 

Now, let drjdt be the radial growth rate of a particle (volatile droplet or soot coated 

dr 

with H2SO4-H2O) in size bin k due to condensation and V. ond = v k + 4k r k — -h be the 

dt 

new particle volume, the rate of change of particle size due to condensation ( V cond >v k ) 
and evaporation ( V cond < v k ) per unit time and per cubic meter of carrier gas can be 
expressed as 

M b 

G k =(£fc,j, k Nj-N k )/h ( 44 ) 

2=1 

where f ckJ , a modified volume splitting operator, is defined as 
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_/+! y cond 

Vj<V C o nd <v j+ 1 

J<M IS 

+ 

1 

v.? 

1 

^ fc,k,j - 1 

V H< V cond< V j 

j> 1 

1 

Vend * Vj 

II 

§ 

ba 

^ 7 cond 

Vcond < V t 

7 = 1 


( 45 ) 


The first term on the right hand side of equation (44) represents the production of 
particles in size bin k due to condensation and evaporation. The last term is the 
destruction of particles in size bin k by condensation and evaporation. 

(dN k ldt) condlevap is expressed as 


dN k ' 


dt 


- G, 


y cond.levap. 


( 46 ) 


The change of radius due to condensation or evaporation is 


_ y 

dt ^ 


dr. 


dn. 


dn.. 


'4 J 


dt 


m- H n O,HnSO, 


where 


dr k _ MW m 
d n m ,k ^Pso! r k 


( 47 ) 


( 48 ) 


The condensation and evaporation processes control the weight fraction of H 2 SO 4 in 
the liquid solution (W H SOi ), on which many thermodynamic quantities depend, e.g., the 
saturation vapor pressure in equation (40). Therefore, it is necessary to determine the 
evolution of W H so . Starting from its definition 


IV, 


MW H 2 SO t * n H 2 S0 4 ,k 


H,SO,,k 


T m MW„ 


x n , 


m = H 2 0,H 2 S0 4 ( 49 ) 


where n m k is the mole number of species m in the droplets of size bin k, and the mass 
of the liquid droplets in the same bin is ^ MW x n m * = 4trp sol r k /3 . Differentiating 
equation (49) with respect to time and invoking equations (47) and (48) leads to 

d^H 2 S0 4 ,k _ 3 
dt r k 


dr. 


\dn H ^ SOi k j 


dn 


H 2 SO..k 


dt 


-IV, 


H.SO,.k 


d _k_ 

dt 


( 50 ) 


Using a semi-implicit formula, the result is 


W, 


IV, 


t-h 

H 2 S0 4 ,k 


H ? SO d ,k 


+ h± - 

r k |_< 

1 + h 


dn. 


h 2 so 4 ,k 


dn 


( 51 ) 
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4.7. Effects of Microphysical processes on Species Evolution 


The conservation equations of gas-phase H 2 SO 4 , SO 3 and H 2 O will have additional 
source/sink terms due to the microphysical processes described above, i.e., 


dC, 


dt 


' WH-yO ' 


^ honft tot 




2 . dn 


H.O.i 


dt 


N 


H.SO.-H.O.i 


T>, 


dn 


H.O.i 


dt 


N soot, i 152) 


dC„ 


dt 


Mb dyi 

. m j * * un H 2 SO A ,i 

Wh 2 S0 4 J hom n tot X 2—1 j, ” H 2 S0 4 -H 2 0.i 


i=l dt 


M B /Jia 

Z A un H 2 S0 4 ,i 

e '^r 


N 

soot,i 


ts 

-.25 acC H , SOi Z ^d sooti N sgotj 


x(l - d t ) 


(53) 


dCs 0 

dt 


' Wso. 


-.25 acC. n V ;r<f 2 f N 

S0 3 / -/ soot, 1 soot,i 


x(l-6».) 


(54) 


It is noted here that, there also exists other models for taking into account the effects of 
particle microphysical processes on the gas species. For example, in reference [20], 
Yelvington used the following relationship [21] to represent the rate of adsorption of 
H 2 SO 4 by soot per aerosol 


B 


H 2 S0 4 \'k 


(h) = 


47tr k D h so 


Ph „ 


4D u c/4 l(.f c h cn v 


'H 2 S0 4 ' 


c,h 2 so 4 v h 2 so 4 


) + fl(r k +A k ) R u T 


(55) 


It is a function of soot particle radius, temperature, and pressure. In addition, D h so ^ , 


H SOi , v H , SOi , and p H so are the diffusivity, molecular accommodation coefficient, 
gas mean speed, and partial pressure of species H2SO4; and \ and r k are the mean free 


path and radius of soot particle in size bin k . 

Similarly, the rate of condensation of a volatile species, a , per wetted aerosol is 
given by [21] as 


C a (r k )[molls] 


4^jP a 

4 D J(P c , a v a ) + rfl(r k + \) 


Pa P a,k 


(56) 


where a is any volatile species, and p°° ak the vapor pressure above a flat surface. 

Using the vapor pressure above a flat surface ignores the Kelvin effect, hence 
underestimates the timescale for condensation. Since the driving force for condensation 
is ( p a - k ), condensation will not occur if the partial pressure is below the vapor 
pressure. Thus, it is very important that the correct saturation vapor pressure is used. 
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5. Applications 

The numerical platform is a NASA in-house code for simulating turbulent combustion 
in propulsion systems. A detailed description of this code can be found in several 
reports [22,23,24]. Its enhancements for post-combustor trace chemistry modeling and 
simulation are reported by Wey and Liu [25]. Under the current effort, we have further 
implemented the microphysical models for aerosol dynamics discussed in the previous 
sections into this analysis tool. In the following, its application to the turbine 
environment as well as to the components of particle sampling system are discussed. 
Several investigations (e.g., Lukachko et al. [3]) have shown the important role of gas 
phase chemistry along the post combustor flow path in the production of sulfate and 
nitrate aerosol precursors. In particular, it has been suggested that multi-dimensional 
analysis needs be applied to the high-pressure turbine section, while one-dimensional 
simulations suffice for the remaining sections of the internal gas path. 

It has also been known that strong modification in the temperature and pressure 
environment experienced by the particles in the sampling system can lead to different 
interpretations of the measurements aimed at characterizing the particle emissions of 
the exhaust. Therefore, there is a need to conduct higher fidelity computations to 
compliment these measurement studies and to assess the effects of sampling 
methodology. 

Finite rate kinetic mechanisms developed for conditions intermediate between 
combustion and atmospheric chemistry have been used in the present work. More 
specifically, a mechanism involving 25 species and 74 reaction steps [26] was used for 
investigating a film cooled turbine blade; later, an improved mechanism involving 29 
species and 73 reaction steps [3] was used for other applications. All simulations 
require the specification of species concentrations along the inflow boundary of the 
computational domain, we have followed the approach outlined in [3] to provide the 
estimation of these species initial conditions. 

5.1. Film Cooled Turbine Blade 

Film cooling injection is widely applied in the thermal design of turbomachinery, as it 
contributes to meet the requirements for reliability and life cycles of modem gas 
turbines operating at high temperature conditions. The film cooling influence on 
post-combustor trace chemistry has been computationally investigated by Wey and Liu 
[27], and a few selected results are cited here. Figure 1 presents the geometry and the 
grid used in the simulation. The coolant mass flow is about five percent of the 
mainstream hot gas mass flow. Five sets of computed results are assembled together to 
show the differences due to various codes as well as the cooling air injection. Data 
reduction has been performed to yield one dimensional distribution along the axis of the 
flow path by averaging the CFD results. The temperature distributions along the axis of 
flow path are depicted in Figure 2. It indicates that the impact of film cooling on the 
temperature is significant. It also suggests that the evaluation of transport properties 
affects the temperature more for the case with cooling air than the case without cooling 
air. NO distributions are shown in Figure 3. The mass fraction for the case with 
cooling air injection decreases significantly as the flow passes through the turbine vane. 
Figure 4 illustrates the SO 3 distributions. Its mass fraction with cooling air injection 
increases more than that without cooling air injection as the flow passes through the 
turbine vane. 
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5.2. First Stage of a High Pressure Turbine 

As part of the NASA/QinetiQ (formerly DERA) emissions testing effort, 
two-dimensional calculations of the evolution of trace species in the first stage of a 
DERA high pressure turbine (HPT1) using the CNEW code were carried out by 
Lukachko et al. [3]. A relatively detailed description of the stage configuration and the 
specification of the inflow condition can be found in [3]. Under the present effort, the 
NCC code was used to calculate the evolution of trace species in the same HPT1 
operating at two nominal power settings - cruise and maximum power - with similar 
high fuel sulfur levels. The results of the high sulfur cruise case will be briefly 
described here; additional results can be found in [28]. Figure 5 outlines the setup of the 
simulations, in particular, the transfer of the solution near the exit of the nozzle guide 
van (NGV) to the inlet of the rotor. Due to the wake of the NGV and the rotation of the 
rotor, the inlet flow of the rotor is non-uniform and unsteady. Figure 6 to Figure 9 
present the fields of the temperature, NO, NO2, and SO3, respectively. In addition to the 
contour plot which is a snapshot at some instant, the area-averaged values also are 
presented. These averaged values are steady in the NGV region but change with time in 
the rotor region. It is noted here that results of the CNEWT code were obtained by using 
the triangular meshes, and only the solutions at one time station are depicted in these 
figures. On the other hand, quadrilateral meshes were used in the code NCC, and 
solutions at several time stations are given. Depending on the variable, these temporal 
solutions may bunch together (e.g., T) or separate more clearly from each other (e.g., 
NO). In general, the solutions obtained by these two codes are quite comparable, except 
for the variable SO3, the reason for this larger difference is not clear. 

5.3. Particulate Matter Probe 

Two particulate matter probes have been investigated: a low dilution probe and a high 
dilution probe. The diluent is N 2 . For the low dilution probe, the tip portion of the inlet 
has an inner diameter of 0.0625 inch, followed by a short divergent section; the 
diameter at the end of the inlet is 0. 1875 inch, the length of the inlet is 0.6553 inch. For 
the high dilution probe, its inlet is a straight tube with an inner diameter of 0.044 inch, 
and a length of 1.2875 inch. In the present work, the computational domain of the low 
dilution probe is set to be 9.558 inch long, and the computational domain of the high 
dilution probe is chosen to be 6.456 inch long. The flow field of the low dilution probe 
was analyzed with an inflow total pressure of 49699 Pa, and inflow total temperature of 
835.78 K. The mass flux of the N 2 diluent is 48 kg nf 2 s~ l at a temperature of 423 K. 
The exit pressure of the probe is set to be 19545 Pa, the wall temperature is set to a fixed 
value of 423 K. The computed dilution factor by mass is about 6.9, which is at the high 
end of the operation of a low dilution probe. Figure 10 indicates that a stagnation zone 
appears in the downstream region of the merged inlet stream and the diluent, while 
another flow stagnation zone sits next to the divergent wall of the inlet tube. These flow 
features will adversely impact the sampling functionality of the probe. Subsequently, 
the same flow conditions were used to evaluate the high dilution probe. However, due 
to the difference in size from the low dilution probe, the computed dilution factor by 
mass now is 25.8. The velocity vector fields of these two probes are compared in Figure 
1 1 . The mixing between the sampled stream and the diluent in these two probes are 
compared in Figure 1 2. The mixing is illustrated by the mass fraction of the sampled 
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stream in the mixture. It is evident that the high dilution probe is much more suitable for 
the given sampling condition. 

Since the high dilution probe has been extensively used in the Aircraft Particle 
Emissions experiment (APEX), additional simulations have been performed for the 
high dilution probe. Conditions for the simulations are listed below: 

• The engine operates at 65% power setting, and the probe is located at 1 
meter downstream of the engine exit. 

• A post combustor trace chemistry kinetic mechanism having 29 species 
and 73 reaction steps is used. 

• Soot entering the probe has a number density of 1.0 x 10" particles per 
cubic meter, it has a log-normal distribution with a median diameter of 40 nm 
and a standard deviation of 1.5. 

• At the entrance of the probe, the pressure of the sampled gas is 101350 
Pa, the static temperature is 720 K, and the speed is 210 meters per second. 
Therefore, the total pressure is 112515 Pa, and the total temperature is 742 K. 

• The static temperature of the diluent N2 is 290 K, the diluent mass flux 
is 1 10 kgm 2 s~\ and the targeted dilution factor by mass is 24. 

• The static pressure of the mixture at the exit of the computational 
domain is 81000 Pa. The temperature of the probe wall is set to 320 K. 

• Aerosol size distribution is divided into 12 bins ( M n = 12 ), the diameter 
of particles in the smallest bin is 3 nm. Soot particles, H 2 SO 4 -H 2 O droplets, and 
the mixture of them are included in the calculations. Consequently, a total of 72 
equations are solved for each mesh cell. 


Since the particle size distribution at the entrance of the probe is assumed to be 
log-normal, the following formula is used to discretize the continuous function into 
stepwise function, for size bin k: 


N 

_ T total 

k Jlirlrn 


-exp 


(lnd pk -lnd pg y 

2 In 2 a „ 
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p,k 
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(57) 


p,k 


where N total is the total number of particles of all sizes per m ' of the exhaust gas, <7, is 
the geometric standard deviation, d p is the median diameter of particles in meter, 
d p k is the diameter of the particles in size bin k . 


The H 2 S0 4 -H 2 0 nucleation rate depends highly on temperature, relative humidity, 
and H 2 S0 4 concentration. The following empirical criterion [29] is used to decide 
whether H 2 S0 4 -H 2 0 nucleation commences: 


C crit - 0.16exp(0. IT - 3. 5RH - 27.7) (58) 

where T is in degree K, RH is the relative humidity on a scale from 0 to 1, and C crit is in 
jug m 3 . When the gas-phase concentration of H 2 S0 4 exceeds C crit , the nucleation of 
H2SO4-H2O droplets starts. Sample results from these simulations are presented in 
Figure 13. While the number density of soot with 66 nm diameter remains essentially 
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unchanged along the probe length, the results indicate the emergence of H2SO4-H2O 
droplets due to nucleation. 


5.4. Sampling Line 

Three-dimensional calculations have been conducted to capture the potentially 
important secondary flow effects by dividing the long sampling line into multiple 
shorter segments. Starting from the exit plane of the high dilution probe, a segment of 
the sampling line is simulated at first; the results at the exit plane of this segment are 
then used as the inlet profiles for the next segment, and so on. Two configurations have 
been investigated. The first one consists of two segments of a simple straight tube. The 
second one consists of a straight tube segment followed by a Z-shaped tube segment. 
Figure 14 presents the temperature distribution along the center line of the probe as well 
as along the center line of the sampling line segments. The temperature inside the 
probe drops rapidly because of the dilution, while the temperature in the sampling line 
is very close to the imposed wall temperature which is kept at 310 K. In Figure 15, the 
mass fraction distribution of species H 2 SO 4 is shown. Due to a sudden expansion in 
volume from the end of the probe to the beginning of the sampling line, a rapid change 
is observed (the typical internal diameter of a sampling line is about 0.34 inch, the 
internal diameter at the end of the currently used AEDC/NASA high dilution probe is 
0.249 inch). In Figure 16, the evolution of the soot particles with 33 nm diameter is 
shown. The loss of soot particles in the straight segment is smaller than the loss in the 
Z-shaped segment. Similar behavior is observed in Figure 17 for soot particles with 66 
nm diameter. Figure 18 suggests that the process of H 2 SO 4 -H 2 O nucleation as indicated 
by the emergence of droplets with 1 nm diameter starts from the inside of the probe and 
continues into the sampling line until the mass fraction of the gaseous H 2 SO 4 levels off 
(see Figure 15). The results in Figure 19 suggest that the coagulation of H 2 SO 4 -H 2 O 
droplets as indicated by the emergence of 2 nm droplets begins in the sampling line. 


5.5. Pressure Reduction Chamber 

The pressure of the exhaust gas at the exit of the combustor is much higher than the 
pressure in the ambient environment. A device which reduces the pressure of the 
sampled exhaust gas is often employed to meet the operating range of the particulate 
matter analyzers. A schematic of a typical pressure reduction vessel is shown in Figure 
20. It is a cylinder-like device that has an inlet tube on the top and an exit tube in the 
bottom. The internal (inside) diameter, ID, of the inlet tube is denoted by a. The portion 
of the inlet tube inside the pressure reducer may be expended up to five degrees to slow 
down the drop of the incoming pressure and the temperature. The gap between the exit 
of the inlet tube and the entrance of the sample extraction tube is denoted by d. The 
diameter of the cylinder is denoted by S. In the current study, the internal diameter of 
the inlet tube is a = 0.09 inch and the diameter at the exit of the inlet tube is .2646 inch 
due to a five-degree expansion. The gap distance, d, is 2 inch; and the ID of the sample 
extraction tube is 0.34 inch. The diameter of the reducer is S = 3 inch; the overall length 
of the domain which includes the inlet tube, the chamber and the sample extraction tube 
is 35 inch (0.889 m). In the current study, instead of conducting three-dimensional 
calculations, axisymmetric simulations have been carried out by replacing the bleeding 
hole with a bleeding slot having the same opening area. Two cases are presented here to 
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show the chemical conversion and particle evolution in the device. The total pressure of 
the incoming gas are 65 and 255 psi, respectively. The total temperature is set to be 485 
K for both cases. The shortest flow residence time is 21 ms for the case of 65 psi and 6.2 
ms for the case of 255 psi. The exit pressure is set at 19.7 psi. The pressure at the 
bleeding location is 14.7 psi. The size distribution of the soot particle is prescribed to be 
log-normal, with median radius of 40 nm, modal widths of 1.5 and total number density 
of 10 13 particle/m 3 . The chemical composition of the inlet gas and the chemical 
kinematic mechanism are the same as those used in the study of particulate matter 
probe. The wall temperature of the inlet tube entering the vessel chamber is set to 400 F 
(477 K). The rest of the wall is either at a temperature of 450 K or insulated. The 
temperature distributions along a line which is 0.018 inch away from the center axis are 
shown in Figure 2 1 . For the 65 psi case, the temperature is above 450 K from the inlet to 
the exit, hence the occurrence of the nucleation of H 2 SO 4 -H 2 O droplets is not possible. 
For the 255 psi case, the temperature has dropped to near 250 K at the end of the inlet 
tube. The distributions of FFSCUmass fraction along the 0.018-inch line are shown in 
Figure 22. The increase of the mass fraction occurs in the straight portion of the inlet 
tube outside of the chamber region. It also indicates that the chamber has little influence 
on the evolution of H 2 SO 4 . In Figure 23, the number density distribution, in logio scale, 
of the soot particles is shown for the 65 psi case. The influence of the chamber on the 
size distribution is very mild. Figure. 24 is for the 255 psi case, the influence of the 
chamber on the size distribution of the soot particles is more noticeable. Figure 25 
shows that, due to the relatively large dropping of the temperature, the nucleation of 
H 2 SO 4 -H 2 O droplets is significant in the case of 255 psi, in addition, the coagulation of 
smaller droplets to form the larger ones also occurs. 

6. Summary 

To better understand the formation and subsequent development of gaseous pollutants, 
particulate pollutants and their precursors of aviation-sourced emissions, it is necessary 
to adopt approaches which combine the use of accurate measurement and high-fidelity 
modeling and simulation. The development of such a CFD based multi-dimensional 
modeling and simulation tool which accounts for and integrates fluid dynamics, 
chemistry and particle microphysics relevant to the aircraft emissions is elucidated in 
the present paper. Its capabilities have been demonstrated by the results of several 
applications in environments representative of post-combustor flow path and sampling 
system. This represents a first step towards the long-term goal of establishing 
physically based tools that can be used to handle fundamental processes acting on 
length/time scales relevant to practical sampling and measurement. Obviously, the tool 
developed under the current effort must be further validated by data from systematic 
and careful parametric measurements. At the present time, measurements aimed at 
providing comprehensive data for the purpose of model development/improvement and 
validation are lacking. This is one area which needs to be strengthened by the research 
community. In addition to validating and further improving the existing predictive 
capabilities, our near term efforts will also include the modeling and simulation in the 
jet engine plume environment and the modeling of soot formation. 
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Appendix A — Parameterized Homogeneous Nucleation Rate 

Because the classical theory is often unreliable for predicting the nucleation rates, 
nucleation parameterizations have been developed to replace the classical theory. 



where In is the nature log function, N sa is the number density of gaseous sulfuric acid 
( molecJcm 3 ), T is the absolute temperature in K , and RH is the relative humidity in 
percent (1 00 xp //2( /p—). 

The nucleation rate is given by an exponential of a third-order polynomial of ln(^-) 
and In (TV J , 


J ho J panicle cm 3 s '] = exp{ a + b - In 




( RH^ 

+e • In (N„ ) + / • In ^ J ln ( N sa ) 

r x m2 

( RH i 

+*• ln [^ jqq J ta(JV„)+A-[ta(JV,)]- 

+M"[^][ ln( A'„ ) ] 2 +y- [ln< A', ) ] ! } 

where the coefficients a...j are functions of temperature and critical cluster mole 
fraction x* 

a = .14309 + 2.21956 • T - .027391 1 • T 2 + .00007228 1 1 • T 3 + 5.91 822/** 
b = . 1 17489 + .462532 T-. Oil 8059 • T 2 + .0000404196 • T 3 + 15.7963/x* 
c = -.215554- .0810269 • T + .00143581 • T 2 -4.77580xl0~ 6 • T 3 - 2.91297/x* 
d = -3.5885600 + .04950800 • T - .0002138200 T 2 +3.10801xl0~ 7 T 3 - .0293333 lx* 
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e = 1. 1459800 - .600796 • T + .00864245 • T 2 - .0000228947 • T 3 - 8.44985 lx* 
f = 2. 1585500 + .0808121 • T - .000407382 • T 2 - 4.01957 x 10~ 7 • T 3 + ,721326/x* 
g = 1.62410 - .01601060 • T + .0000377124 • T 2 + 3.21794 x 10 -8 • T 3 - .01 13255 lx* 
h = 9.7 16820 - . 1 1 50480 • T + .000157098 • T 2 + 4.009 14 x 10^ 7 • T 3 + .7 1 1 86/x* 
i = -1 .0561 1 + .00903378 • T - .0000198417 • T 2 + 2.46048 x 10 -8 • T 3 - ,0579087/x* 
j = 148712 + .00283508 • T - 9.246 19 d -6 T 2 +5.00427 xl0“ 9 T 3 - .012708 1 lx* 

The total number of molecules in the critical cluster, N * ot , is given by 

r f RH ^ 

N, n . = expl ca + cb - In 

1 UooJ 



r ( rh Y 
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r ( rr \] 
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+ cd ■ 
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UooJ 
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+ce ■ ln( N sa ) + cf Tn 


RH^j 

Too, 


1 n (AU 


+cg- 


In 


RH 


100 


\n(NJ + ch-[MNj] 2 


+ci ■ In 


RH 

loo 


HNj] 2 + cj-[ln(Njf } 


where the coefficients ca...cj are functions of temperature and critical cluster mole 


fraction x* 

ca = -.00295413 - .097683400 • T + .00 10248500 • T 2 - 2. 1 8646 x 1 O’ 6 • T 3 - . 10 17 17000/x* 
cb = -.00205064 - .007585040 • T + .0001926540 • T 2 - 6.70430 xl0“ 7 T 3 - ,255774000/x* 
cc = 3.22308 d - 3 + .000852637 T - 1 5A151d -5 T 2 + 5.66661 x 10 -8 • T 3 + ,03384440/x* 
cd = .047432300 - .000625104 • T + 2.65066 d -6 T 2 - 3.6747 lxlO -9 T 3 - .00026725 1/x* 
ce = -.0125211000 + .005806550 • T - .0001016740 T 2 +2.88 195 xlO” 7 -T 3 + ,094224300/x* 
cf = -.0385460000 -.0006723 16 T + 2.602SSd -6-T 2 + 1. 19416 x 10“ 8 -T 3 - ,008515150/x* 
eg = -.018374900 + .0001720720 • T - 3.71766rf -I T 2 - 5. 14875 x 10 10 • T 3 + .000268660 lx* 
ch = -.061997400 + 9.06958 d -AT -9.1 1728d -IT 2 - 5.36796 x 10 -9 • T 3 - ,007742340/x* 
ci = .012182700 - .0001066500 • T + 2.53460r/ - 7 • T 2 - 3.635 19 x 10 10 • T 3 + .00061 0065/x* 
cj = .000320184 - .0000174762 • T + 6.06504t/ - 8 • T 2 - 1.42177 x 10" 11 • T 3 + .00013575 1/x* 


The radius of the cluster in nanometers is given as a function of the mole fraction and 
the total number of molecules in the cluster: 


r [ nm ] = exp [- 1 .6524245 + 0.423 1 6402x* + .33466483 MNf )] 

A.2. Parameterization of High Temperature Nucleation 

This parameterization is valid for temperature ranging from 300.15-400.15/^ and 
relative humidity ranging from 1-100% . The mole fraction of sulfuric acid (or the 
acidity) in the critical cluster is given by 
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x*=. 847012-0.00296567 -0.00662266 ln(/VJ + .00005878357Tn(/VJ 
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where N sa is the number density of gaseous sulfuric acid ( molecJcm 3 ), T is the 
absolute temperature in K , and RH is the relative humidity in percent 
( 1 00X PhjJPh.,0 )• 

The nucleation rate is given by an exponential of a third-order polynomial of ln(-SS-) 
and In (N J 


J honXP art ^ e cm ' s ' I = exp{ a + b In 


RH^| 

Too, 



r ( rhY 
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\n(NJ + h-[\n(Nj] 2 


HNj] 2 + j-[ln(Njf } 


where the coefficients a...i are functions of temperature and critical cluster mole 
fraction x* 


a — -.00156975 - .134245- T + .100507 • T 2 -.000460 103 -T 3 + .18741 6/ x * 2 + . 01041 22/x* 
b = .00195077 + . 168038 • T - .0225755 • T 2 + .0000827149 • T 3 + .0025029/ x * 2 + ,0155215/x* 
c = .00015408 - .0280301 • T + .00154587 • T 2 - 4.5270 lxlO" 6 T 3 + .091 5323/ x * 2 + .07 1 1 652/ x* 
d = -.00509267 - .00796846 • T + .0000446828 • T 2 - 8.79425 x 10 -8 • T 3 + . 133991/ x * 2 + .831 1 12 /jc* 
e = -.0227223- 1.56512-7’ + .00380717 -T 2 + 1.64109//- 5- T 3 + 1.29499/ x * 2 + .047482 1/x* 

/ = . 00310646 + .304518 T - . 000564012 -T 2 -2.03267xl0" 6 -T 3 -,351584/ x * 2 + .103749/x* 
g = .077543-. 00196315- T-. 0000130412 T 2 + 6.62369xl0" 8 -T 3 +. 01 1347/ x * 2 + .0972804/x* 
h = -.l 53 143 + .0575392 • T - .0003065 1 1 • T 2 - 2.96097 x 10 -8 • T 3 - .09825 14/ x * 2 + ,336286/x* 
i = -.552173 - .00207043 • T + .0000144032 • T 2 + 8.83 x 10“ 9 • T 3 + .01 19833/ x * 2 - ,0700025/x* 
j = .126544 - .00136029 • T + 5.90598 d - 6 ■ T 2 - 4. 17 15 x 10" 9 • T 3 + ,00170807/ x * 2 - ,0064323/x* 


The total number of molecules in the critical cluster, N * ot , is given by 
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where the coefficients ca...cj are functions of temperature and critical cluster mole 


fraction x* 

ca = 7.5 1024 d - 6 + .000502054 • T - .0000368602 • T 2 + 1 .08256 x 10" 6 • T 3 - ,000270282/x* 
cb = -4.30048 d - 6 - .000730133 • T + .0002520620 • T 2 - 1 .01648 x 10“ 6 • T 3 - .001 142830/x* 
cc = -4.42 156 d - 6 - .002348600 • T + 3.0065 xl0“ 7 - T 2 + 2.44797 xl0“ 8 -T 3 - .002502260 lx* 
cd = -.000167057 + . 000207504- T -1.13013xl0" 6 -T 2 + 1.80268xl0 -9 -T 3 - .016824500/x* 
ce = .0000985954 + .004512850 • T - .0000512557 T 2 + 4.60749 xl0“ 8 -T 3 - .002143 1 80/x* 
cf = .0000636528 - .002885290 • T + 6.51706 x 10 -6 • T 2 + 2.32601 xl0“ 8 • T 3 - .001 103 190/x* 
eg = .000449239 + .0000689416 • T - 3.50302 x 10' 7 • T 2 + 1 .0745 1 x 10 MO • T 3 + .001696460 lx* 
ch = .00083 1 844 - 5.35 108 x 10" 6 • T + 1 .66432 x 10" 6 • T 2 - 3.05 108xl0“ 9 T 3 - .00030625 1 lx* 
ci = .003553744 + .0000306009 • T - 2. 1 1004 x 10” 7 • T 2 - 2. 1 1436 x 10~ n • T 3 + ,000749890/x* 
cj = -.00 1435340 + 7. 85600 x 10“ 6 - T - 3.45 128 xlO -8 -T 2 + 5. 21547x10^' -T 3 - .00002 1423/x* 


The radius of the cluster in nanometers is given as a function of the mole fraction and 
the total number of molecules in the cluster: 


r[nm\ = exp [-1. 6525507 + 0.45852848x +. 33483673 ln(AT,)] 
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Appendix B — Thermodynamic Properties 

In order to solve the aerosol dynamics, various thermodynamic properties are needed 
for the calculation of the microphysical processes. 

B.l. Saturation Pure Water Vapor Pressure 

According to Preining et at. [30], the saturation vapor pressure ( Nlnr,T in K ) for pure 
water is 


/Co = exp^ 77.34491296- 


7235.424651 


-8.21nT + 0.00571 137A 


B.2. Saturation Pure Sulfuric Acid Vapor Pressure 


According to Kulmala and Laaksonen [31], the saturation vapor pressure 
( N I m 2 ,T in K ) for pure sulfuric acid above a flat surface of H 2 SO 4 liquid is 


v flat 

y0,H 2 SO< 


= 101325.0xexp{ — 11.695 + 10156xf 

T 360.15 


+ 


0.38 

545" 


x l. + ln( 


360.15 ) 360.15 


] } 


The range of the temperature is from 190. 15 K to 400. 15 K. 


B.3. Surface Tension of Sulfuric Acid and Water 

The following fit gives the surface tension o(J ! nr ) for all sulfuric acid solution as a 
function of the temperature T and the mole fraction x : 

<j(x,T)\Jm~ 2 ] = (a + bT l )- 7 ] 1 256 

where 


a{ x) = .2358 - .529 • x + 4.073 • x 2 - 12.6707 • x 3 +15.3552 • x 4 - 6.3 138 • x 5 
b(x) = -.14738 + 0.6253 • x - 5.4808 • x 2 + 17.2366 • x 3 - 21.0487 • x 4 + 8.719 • x 5 

Tfx,T) = \-j 

The pseudo critical temperature T c (x) is given by 

T c (x) = 647.15- (1.0- x ) 2 + 900.0 • x 2 + 3 1 56. 1 860 • x • (1 .0 - x) 

Note that the extrapolation to higher or lower temperatures behaves smoothly. 

B.4. Density of Sulfuric Acid Solution 

For temperature ranging from 273 K to 373 K, the density p sol (kglm 3 ) of the sulfuric 
acid solution is fitted to a polynomial function of temperature and sulfuric acid mass 
fraction x m [32]: 

P sol (x m ,T) = ( a(x m ) + b(x m ) ■ T + c(xj • T 2 )- 1000.0 
Here T is the temperature of the liquid ( K ) and x m is the mass fraction of sulfuric acid. 
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The coefficients a(x m ) , b( x m ) and c(x m ) are obtained from the following equations: 
a(x ) = .7 68 1724 + 2.1 847 140 • x +7.1 63002 • x 2 m - 44.3 1 447 • x 3 
+88.75606-x 4 -75.73729-x„ 5 ,+23.43228-x„ 6 , 


b( x ) = 1.808225 • 10" 3 -9.294656 • 10 -3 ■ x - .03742148 • x; 
+0.2565321- x 3 -.5362872- x 4 +.4857736- x 5 

m m m 

-.1629592-x 6 

m 

c(x ) = -3.478524 • 1CT 6 + 1.335867 • 1CT 5 • x„, 

+5. 195706-10~ 5 -x,;-3. 717636-1 CH-x 3 
+7.9908 1 1-1 CH-x 4 -7.45 8060-1 0" 4 -x 5 

m m 

+2.58139-10 _4 -x 6 

m 


B.5. Partial Molar Volume in a Two-Component Solution 


The Kelvin equation involves v„, ( nr’hnol ) , i.e., the partial molar volume of species 
m in a two-component solution. Given the volume v and the mole fraction of 
component 2, i.e., x 2 , the partial molar volumes pi and y 2 are computed from the 
volume change when mixing the pure components: 

_ AV 

Vi = V x 2 

AX 2 


and 


V 2 = v (x 2 — 1) 

AX 2 

such that (1 -x 7 )v 1 + Xtv 2 = v . The molar volume of a two-component solution is 
given by 

- _ MW s„, _ (1 ~ vWWi + VMW 2 (B1 

P:-f[ P sol 

The slope, , is computed by using the finite difference: 

av V-v 


A^2 X 2 2 

where v" and v' are evaluated by slightly perturbing x, in Equation (Bl). 


B.6. Saturated Partial Pressure of H 2 S0 4 and H 2 0 Vapors 

The Kelvin equation also involves P n f ’ a s ' at , i.e., the saturated vapor pressure of species m 
over a flat liquid surface having the same composition as the droplet. It is defined as 

p!!“L = r,n x m p 0 , m > m = h 2 o,h 2 so a 


where y m is the activity coefficient, p 0m is the saturation vapor pressure of pure 
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component m . The activity coefficients, y HO and y H vo _ , in a binary system can be 
expressed as 


t- log r Hl0 


^H,0 X 


2 

h 2 so 4 


[x t 


h 2 so 4 ^h 2 o x h 1 o^ 


T ■ log y, 


A Y Z 

ri H 2 S0 4 A 'H 2 0 


H 2 S0 4 r . n -i2 

[ X //,0 + B H,S0 4 X H 2 S0 . ] 


where T is the temperature in degree K , x is the mole fraction, and the coefficients 
A and B are given by Zeleznik [33]: 


A HiQ = 2 . 989 x 10 


3 2.147x10° 2.33x10 


- + - 


4 ,n 3 4.074x10° 4.421x10 s 

A a ™ =5.672x10 


1 h 2 so 4 


- + - 


Bh 2 o 


B 


■ = 0.527 


// 2 S0 4 


5.7. Diffusivity of H 2 S0 4 and H 2 0 Molecules 


According to Pruppacher and Klett [34], the diffusivity of H 2 O molecules ( m 2 ls ) at 
temperature between 233 K and 313 K is given by 


D 


v,H 2 0 


= 10x0.211 — 




1 7 noA 


.94 




where T" = 273. 1 5K , P" = 1013.25 mbar. Since D v <xvA g and the mean thermal 
velocity v °< 1 1 yj MW v > oc 1 1 yj MW v ■ Therefore, D vM so can be obtained through 

JmWh 2 o 


P^v,H 2 SO. x 


Jmw 


h 2 so 4 


Alternatively [35], d, h so can be computed from 


D 


8 Ad 2 p 


-x 


RJ MW, 
2k 


MW h 2 so 4 + MW g 
MW h 2 so , 


where d is the diameter of a sulfuric acid molecule, commonly taken as 4.5x10 10m , 
MW h, so, ' s 4s molecular weight ( kg/mole ), MW g ' s the mass of carrier gas per mole 

( kglmole ), p is the density of carrier gas (kg Ind) and A is the Avogadro’s number 

(= 6.02252 xl0 23 mo/c _1 ). 
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A TURBINE VANE WITH 12 COOLING SLOTS 



Figure 1: A two-block overset grid around the vane. (Thickness of vane is about .01m.) 


NASA/TM- 


- 2006-214397 


30 


Film Cooling Effect on Temperature for A Generic Turbine Vane 


Area- Averaged Temperature Distribution Along Axis of Flow Path 



Figure 2: Comparison of area-averaged temperature distributions along the axis of flow path. 
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NO (MASS FRACTION) 


Film Cooling Effect on NO for A Generic Turbine Vane 

Area- Averaged NO Distribution Along Axis of Flow Path 



Figure 3: Comparison of area-averaged NO distributions along the axis of flow path. 
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S03 (MASS FRACTION) 


Fi lm Cooling Effect on S03 for A Generic Turbine Vane 


Area- Averaged S03 Distribution Along Axis of Flow Path 



Figure 4: Comparison of area-averaged S0 3 distributions along the axis of flow path. 
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Figure 5 Boundary condition specification for 2-D unsteady HPT1 simulations 
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NGV and Rolor of DERA High Sulfur Cruise Case 




Figure 6 Area-averaged temperature distribution along the axis of flow path 
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NGV and Rolur of DERA High Sulfur Cruise Case 
Area-Averaged NO Distribution Along Axis of Flow Path 
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Figure 7 Area-averaged NO distribution along the axis of flow path 
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N02 (MASS FRACTION) 


NGV and Rotor of DERA High Sulfur Cruise Case 

Area-Averaged N02 Distribution Along Axis of Flow Path 
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Figure 8 Area-averaged N02 distribution along the axis of flow path 
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S03 (MASS FRACTION) 


NGV and Rotor of DERA High Sulfur Cruise Case 
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Figure 9 Area-averaged S03 distribution along the axis of flow path 
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Stream 


NASA/AEDC low dilution PM probe 

Stream: P tot ai=49699 N/M**2, T to t a i=835.78 K; 

Exit: P eX it=19545 N/M**2; Diluent N 2 : Speed=150 M/S, 
density .32 KG/M**3, T sta tic =423 K; T wa „=423 K 
Dilution factor by mass =6.15 



Streamlines are colored by Mach number 


Figure 10 Streamlines of the low dilution PM probe 
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Stream: P_tot=49699 N/M**2, T_tot=835.78 K; Exit: P_static= 19545 N/M**2 
Diluent N2: Speed=150 M/S, Density .32 KG/M*’3, T_static=423; 

T_wall=423 K 


NASA/AEDC Low Dilution PM Probe 
^.Calculated Mass Ratio = Diluent N2/ Stream Air = 6.9 






VELOCITY VECTOR (M/S) 

0.02642 161.9 323 6 485 8 647.7 809.6 



NASA/AEDC High Dilution PM Probe 

Calculated Mass Ratio = Diluent N 2 / Stream Air = 25.8 



VELOCITY VECTOR (M/S) 

0.003915 65.27 130.5 195.8 261.1 326.3 



Figure 11 Velocity vectors of two different PM probes 
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Stream: P_tot=49699 N/M**2, T_tot=835.78 K; Exit: P_static= 19545 N/M**2 
Diluent N 2 : Speed=150 M/S, Density .32 KG/M*’3, T_static=423; 

T_wall=423 K 


NASA/AEDC Low Dilution PM Probe 
Calculated Mass Ratio = Diluent N2/ Stream Air = 6.9 



Stream Mass Fraction 


0 0.2 0.4 0.6 0.8 1 



NASA/AEDC High Dilution PM Probe 

Calculated Mass Ratio = Diluent N2/ Stream Air = 25.8 



Stream Mass Fraction 



I 


Figure 12 Mixing between the sampled stream and the diluent. 
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Number Density of Particles in High Dilution Probe 

Pt=112515 pa, Tt=742 K, P_exit=81000 pa, T_waU=320 K 



Figure 13 Soot evolution and H2S04-H20 nucleation in high dilution probe. 
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High Dilution Probe with Sampling Line 

Temperature Distribution Along the Center line 


n — i — i — r~ 


t 

V 


“i — i — i — r~ 


n — i — i r~ 


n — i — i — r 


n — i — i — r 


G-© High Dilution Probe 
□-0 First Segment (Straight) 
o-€> Second Segment (Straight) 
Second Segment (Z-shape) 


E- DDDatuaDt nmuji tt 


J I I I I I I I I I I I I I I I I I I L 


bUr fi 


0 -8 -6 -4 -2 0 

Center Line of Probe and Sampling Line (meter) 


Figure 14 The temperature distribution along the center line of the probe and sampling line. 
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High Dilution Probe with Sampling Line 



Center Line of Probe and Sampling Line (meter) 


Figure 15 The gaseous H 2 S0 4 distribution for the probe and sampling line. 
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High Dilution Probe with Sampling Line 

Number Density of Soot (diameter ~ 33 nm) along the Center Line 



Center Line of Probe and Sampling Line (meter) 


Figure 16 The number density distribution of soot particle (33 nm) for probe and sampling line. 
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High Dilution Probe with Sampling Line 



Center Line of Probe and Sampling Line (meter) 


Figure 17 The number density distribution of soot particle ( 66 nm) for probe and sampling line. 
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High Dilution Probe with Sampling Line 



Center Line of Probe and Sampling Line (meter) 


Figure 18 The number density distribution of H 2 S0 4 -H 2 0 droplet for probe and sampling line. 
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High Dilution Probe with Sampling Line 



Center Line of Probe and Sampling Line (meter) 


Figure 19 The evolution of H2SO4-H2O droplet in probe and sampling line. 
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Pressure Reduction Vessel 


Excess Flow 
0.5” OD 


a=Inlet tube 1D=0.09 to 0.125” Machined 



Figure 20 Schematic of a pressure reduction vessel. 
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Pressure Reduction Chamber Study (400 F case) 

Temperature Distribution Along .018 Inch Above of Axis 



Figure 21 The temperature distribution along a line close to the center line of the device. 
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H2S04 (Mass Fraction) 



Representative Flow Path (meter) 


Figure 22 The H 2 S0 4 distribution along a line close to the center line of the device. 
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Pressure Reduction Chamber Study (400 F case) 



Representative Flow Path (meter) 


Figure 23 The number density distribution of soot particles. 
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Pressure Reduction Chamber Study (400 F case) 



Representative Flow Path (meter) 


Figure 24 The number density distribution of soot particles. 
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Pressure Reduction Chamber Study (400 F case) 



Representative Flow Path (meter) 


Figure 25 The number density distribution of H 2 SO 4 -H 2 O droplets. 
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